Cold quarks in medium: an equation of state 
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We derive a compact, semi-algebraic expression for the cold quark matter equation of state (EoS) 
in a covariant model that exhibits coincident deconfinement and chiral symmetry restoring transi- 
tions in-medium. Along the way we obtain algebraic expressions for: the number- and scalar- density 
distributions in both the confining Nambu and deconfined Wigner phases; and the vacuum-pressure 
difference between these phases, which defines a bag constant. The confining interaction materially 
alters the distribution functions from those of a Fermi gas and consequently has a significant impact 
on the model's thermodynamic properties, which is apparent in the EoS. 
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A reliable equation of state (EoS) for cold quark mat- 
ter would be extremely valuable in modern astrophysics. 
As emphasised, e.g., by Ref. [1], the identification of a 
neutron star with a quark matter core depends upon it. 
The problem is broader. An almost complete absence of 
experimental constraints at densities above nuclear sat- 
uration entails that little is truly known about the EoS 
of any form of dense matter beyond that point [2] . 

In attempting to predict the properties and astrophys- 
ical signals of cold quark matter, all we have currently 
are models. Owing, amongst other things, to the so- 
called fcrmion sign problem, the numerical simulation of 
lattice-QCD will not supply this need in the foreseeable 
future [3] . Two classes of models are widely used in this 
application: bag- and Nambu- Jona-Lasinio-models. As 
usually formulated, bag-like models possess a form of con- 
finement but cannot describe dynamical chiral symme- 
try breaking (DCSB), whilst Nambu- Jona-Lasinio (NJL) 
models express DCSB but not confinement. (NB. The 
static potential measured in quenched lattice-QCD is not 
related in any known way to the question of light-quark 
confinement, which can be connected with the analytic 
properties of QCD's Schwinger functions [4].) 

Dyson-Schwinger equations (DSEs) provide a contin- 
uum approach to QCD that can simultaneously address 
both confinement and DCSB. They have been applied 
with success to hadron physics in- vacuum; e.g., [5, 6, 7], 
and to QCD at nonzero chemical-potential and tempera- 
ture [8] . In the context of cold quark matter, the stability 
and properties of various Cooper-paired phases was re- 
cently explored via a truncation of the gap equation [9] 
but that study did not provide information on the EoS. 

An attempt to calculate the EoS was made in Ref. [10]. 
However, in adopting a mcromorphic model for the quark 
propagator and failing to consider the role played by 
the pressure difference between the Wigner and Nambu 
phases, this study could not describe a first-order chiral 
symmetry restoring transition. As exhibited in Ref. [11] 



and elucidated in Ref. [12], coincident, first-order chiral 
and deconfinement transitions are the natural results. 

We will develop an EoS for a Poincare covariant model 
[13], whose gap equation's solution does not admit a 
meromorphic parametrisation, and which exhibits coin- 
cident, first-order deconfinement and chiral symmetry 
restoring transitions. A numerical analysis that exposes 
aspects of this model's EoS is described in Ref. [14]. 
However, we obtain analytic results. In particular, 
for the quark number- and scalar-density distributions, 
fi(\p\ j (J>) an d /2dpi ; At), which will themselves be useful, 
e.g., in astrophysics explorations. Furthermore, these for- 
mulae arc helpful in elucidating novel possibilities for the 
behaviour of fermions subjected to an interaction that 
supports a confining Nambu- Goldstonc phase. 

The in-medium, dressed-quark propagator is [15] 



S(p;n) 1 = ■ p A(p 2 , p ■ u) 

+ i-n{p i + iy)C(p 2 ,p-u) + B(p 2 ,p-u) 



(1) 



where u = (0, i/x), with fj, the quark chemical potential, 
and, in our Euclidean metric: {7 P ,7cr} = 25 pcr ; 7^ = nf p . 
(NB. We employ an ultraviolet-finite model and hence no 
discussion of rcgularisation or renormalisation is neces- 
sary.) The propagator is obtained from the gap equation 

Sip-^y 1 = ij-p + ijiipi +ifi) +771 + £(p;ju) , (2) 



(3) 



where m is the bare mass, D pcr (k;fi) is the dressed- 
gluon propagator and (q, p; /j) is the dressed-quark- 
gluon vertex. 

We specify the model through the choices [13] 



g 2 D pa {k) = S pa 



kp k a 



47rV<5 4 (fc) , 



(4) 
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with 77 a mass-scale parameter, and T®(q,p) = -X a ^/ lT . 

This vertex defines a rainbow gap equation, which is 
the leading-order in a systematic, symmetry- preserving 
DSE truncation scheme [16, 17]. The infrared enhance- 
ment exhibited by Eq. (4) provides for confinement and 
DCSB [5], and the model is super-asymptotically-frcc 
because the interaction strength vanishes for nonzero 
relative momentum. In practice, the model has many 
features in common with a class of renormalisation- 
group-improved effective-interactions; and its distinctive 
momentum-dependence works to advantage in reducing 
integral- to algebraic-equations that preserve the charac- 
ter of the original. It has been used widely with success; 
e.g., in exploring the impact of dressing the quark-gluon 
vertex [18, 19, 20] and in illuminating general, exact re- 
sults connected with the Ua{1) anomaly [21]. 

In the chiral limit, the nonperturbative, chiral sym- 
metry preserving solution of the model gap equation is 
A(p 2 ,p-u) = C(p 2 ,p-u), 



C{p\p-u) 



1 



1 



2'/ 



, B(p 2 ,p-u)=0, (5) 



where p 2 = p 2 + 2p ■ u + u 2 . It describes a phase in which 
chiral symmetry is realised in the Wigner-Weyl mode and 
the quark is not confined. For the analysis which follows, 
it is important to note that C(p 2 ,p-u) possesses a branch 
point at 2if + p 2 + p\ — i-i 2 = 2p 4j tt = 0. For p ^ 0, it 
occurs at p4 = 0, p 2 + 2^ 2 = p 2 . Hence, the branch point 
plays a role when p 2 < 2r] 2 . 

The gap equation also has a confining solution, in 
which chiral symmetry is dynamically broken; viz., for 
m = 0, A(p 2 ,p ■ u) = C(p 2 ,p ■ u), 
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FIG. 1: Wigner-phase quark number density distribution, 
Eqs. (10), (13). Upper panel - p, = 0.5, 1, 2, 4, 8, 16 GeV: the 
Fermi gas result is recovered with increasing p. Lower panel - 
{Q\P) increases with increasing p for p 2 < 2n 2 , illustrated 
at equidistant p G [0.2, 2.0] GeV. At fixed p 2 > 2rj 2 , fi drops 
from 1 to on the domain fi 2 — 2n 2 < jf < /i 2 . 



C(p 2 ,p-u) = 



B(p 2 ,p- 



1 



2 T 



Re(p 2 ) < \ 
otherwise, 



vV 





Ap 2 



otherwise. 



(G) 



(T) 



It describes a phase in which chiral symmetry is realised 
in the Nambu-Goldstonc mode. Confinement is signalled 
by a square-root branch point at p 2 = ?7 2 /4, associated 
with the scalar piece of the self energy. For /1 7^ 0, it 
occurs at pi = 0, p 2 = /i 2 + ?/ 2 /4. 

Equations (5) and (6), (7) are all one needs in order to 
obtain this model's zero-temperature EoS. NB. We ignore 
quark Cooper pairing herein. Diquark condensates have 
previously been considered within this model [22]. We 
will extend our analysis to that case in the future. 

We express the single-quark number density 



n q {n) 



d 3 p 



fi 



(8) 



(2tt) 3 
1 f°° 

— J dp 4 tr D [-745'(p; Mj] , (9) 



where the trace is over spinor indices alone. In the 



Wigncr phase, using Eqs. (5), one obtains 

p 2 < [i 2 - 2?7 2 
p? - 2r, 2 < p 2 < p? (10) 



/r(b1;^) 



i, 



, p 2 < p 2 



where, with A = [p 2 - \p\ 2 ]/[2t] 2 }, 



i w (\p\;p) = 1 + - v/A 2 (l- A 2 )-arccosA 



(11) 



The noninteracting Fermi gas result is recovered from 
Eq. (10) when p ^ rp viz. [upper panel of Fig. 1], 



■^eu-p). 



(12) 



Novel features of the single-quark number density dis- 
tribution are only exposed when p < 77. They origi- 
nate in the dressed-quark propagator's branch point at 
p 2 = p 2 - 2rj 2 [sec Eq. (5)]. 

Whether this curious behaviour is important in the 
consideration of compact stars and heavy ion collisions 
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depends on the natural scale for rj. That scale is set via 
the meson spectrum, a procedure which yields [13] 



i] w ^2m p = 1.09 GeV. 



(13) 



This being so, the unconventional behaviour is important 
and one cannot justify the treatment of quark matter 
as a quasi-ideal Fermi gas; e.g., as a system in which 
the medium serves only to produce a density-dependent 
shift in the chemical potential. Emphasis is found in the 
following facts. For < /u < ^2r\, the noninteracting 
Fermi gas result fi^{\p\]fi) = 1 is precluded. Indeed, 
it follows from Eq. (11) and is illustrated in the lower 
panel of Fig. 1, that on this /i-domain the number density 
distribution is < 1, V|p|. Furthermore, 



n 



0;/j) 



2^2 fx 

7T 7] ' 



(14) 



i.e., at zero momentum the number density increases lin- 
early with chemical potential. 

The Nambu-phase single-quark number density is ob- 
tained from Eqs. (6), (7), (9); viz., 



47r/f(b1;/i) 



dp 4 tr D [-74«S/v(IpI ,P4;aO] 



i->i 



2Re / dp 4 tr D [-74^(1^ ,p4;/i)], (15) 
'pI 



with pi = Re^/j 2 + 7] 2 /4-p 2 . Evidently, only the 2 nd 
term on the rhs contributes for p 2 > /i 2 + 7/ 2 /4 and it 
evaluates to the last line of Eq. (10), so that 



p 2 >li 2 +ri 2 /4 



0. 



(16) 



The Nambu-phase number density evolves in a new 
and unusual manner on p 2 < fi 2 + ij 2 /4; viz., 



4M 




(17) 



The 1 st term is 



where p* 2 — \ + 2i/iy ^ — p' , 
always negative and greater in magnitude than the posi- 
tive, 2 nd term. Hence, fi{\p\ ; ft) < for p 2 < fi 2 + ?? 2 /4 
[see Fig. 2]. Moreover, f^(\ji\;fi ->■ 0+) ->■ 0" and 
f^(\p\;fi — > oo) — > — oo. Notably, even for /i ~ 0, 
the Nambu-phase number density is negative-definite for 
|p] < 77/2. It is plain that the curious features of 
are again closely connected with the branch point in the 
quark propagator [see Eqs. (6), (7)], which in this phase 
is intimately associated with confinement and DCSB. 
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FIG. 2: Nambu-phase quark number density distribution, 
Eq. (17); top to bottom: ft = 0.2,0.4,0.6,0.8,1.0 GeV. The 
density decreases with increasing /1, as it must: the confined, 
DCSB phase is destabilised by increasing /x and hence the 
pressure associated with this vacuum phase must decrease. 



The chiral-limit condensate can be obtained from the 
single-quark scalar density distribution 



1 r°° 

h{m;v) = -rl dp 4 tr D {S(\p\ 7 p 4 ;fj,)} 



(18) 



Since the chiral-limit drcsscd-quark propagator is trace- 
less in the Wigner phase, this density is identically zero. 
On the other hand, in the Nambu phase one obtains 

/f(|pl ; fi) = T h (pi, \p\ ; fi) ~ X h (0, \p\;fi), (19) 

with (z^ = z + ifi) 



4 T] 2 




p 2 + z 2 



— 



(20) 



which is positive-semi-definite, as illustrated in Fig. 3. 
This expression yields 

which is independent of rj, so long as 7/ ^ 0. Since 
the scalar density distribution is dimensionless, an 77- 
independent result had to be obtained in this limit. Nat- 
urally, the scalar density distribution vanishes for 77 = 0; 
i.e., in the absence of interactions. 
The chiral-limit quark condensate is 



(qq)(fx) = 2N c 



* f?(\p\w) 



(27T) 



(22) 



and it is plain from Fig. 3 that the condensate must in- 
crease in magnitude with increasing ft. Moreover, a com- 
parison of Figs. 2 and 3 shows that the evolution of the 
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FIG. 3: Nambu-phase single-quark scalar density distribu- 
tion, f2(\p\ ; n) in Eq. (19). Dashed curve - fi = 0; and solid 
curves - fj. = 0.2, 0.4, 0.6, 0.8, 1.0 GeV. 

Nambu-phase scalar density distribution is anticorrclated 
with that of the number distribution. One may now un- 
derstand that the decrease in the Nambu-phase pressure, 
which can be inferred from Fig. 2, is connected with the 
energy cost of rearranging the vacuum so as to increase 
the magnitude of the condensate in the face of opposi- 
tion from the rising chemical potential. The behaviour 
of hadron properties under these and similar conditions 
is illustrated in Rcfs. [23, 24]. 

Hitherto, general qualitative features of the pressure 
were used in explaining aspects of our results. We now 
present this model's EoS. As the model lies within the 
class of rainbow-truncations, the vacuum-pressure asso- 
ciated with a given phase can be calculated using the 
"steepest-descent" approximation; namely, 

P[S] = TrLn [S" 1 ] - ^Tr [E S] . (23) 

Equation (23) is the auxiliary field effective action eval- 
uated at its minimum [25]. Owing to Eq. (4), in this 
analysis we can neglect the gluon contribution. 

The pressure difference B(ji) = P[S N ] ~P[S W ] is com- 
puted using Eqs. (5), (6), (7). It can be identified with a 
bag constant [26] and, for two light flavours, yields [14] 

B Nf=2 = (0.102 r/) 4 = (0.111 GeV) 4 =: e 4 . (24) 

As evident in Fig. 4, a first-order transition from the 
Nambu- vacuum (confined, with DCSB), to the Wigner- 
vacuum (deconfined, chirally symmetric) occurs at 

/i cr := {(J,\B(jj.) = 0} = 0.276 rj = 300 McV. (25) 

The bag constant receives no contribution from s- 
quarks because their current-mass is too large to support 
a Wigner phase [27]. This does not preclude a Cooper- 
paired phase involving s-quarks at very large /j,. Also, e v 
and ~ 40% larger in models with an interaction 

whose ultraviolet behaviour is more realistic [11, 12]. 

The thermodynamic pressure of each phase is deter- 
mined via the quark number density. We compute it from 
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FIG. 4: Dotted curve - Difference between Nambu-Goldstone 
and Wigner- Weyl vacuum-pressures, which vanishes at /i cr = 
300 MeV; dashed curve - thermodynamic pressure in the 
Nambu phase; Solid curve - thermodynamic pressure in the 
Wigner phase, which is the ground state for /i > /i cr . 
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FIG. 5: Solid curve - Thermodynamic pressure associated 
with the Wigner- Weyl phase of our model; and dashed curve 
- pressure of a noninteracting Fermi gas. Evidently, P? oc fi 5 
for fi < 77 and the chemical potential overwhelms the interac- 
tion for ji > 4 77. 



the generating functional in steepest-descent approxima- 
tion; viz., using Eqs. (10), (16), (17) (<j> = W,N): 

= P^+J\\zni{z),nt^) = ^£dzz 2 f»(z), 

(26) 

where P% N = B - P$ and P% w = -P£, with P* defined 
so that P^(Hc) = 0. P™' N (v) are depicted in Fig. 4. 

Now recall Eq. (14): f^(\p\ = 0) cx \i for i? < 2?/ 2 . 
It follows, using Eqs. (26), that P^ oc /i 5 on this do- 
main. This behaviour stands in marked contrast to that 
of a noninteracting Fermi gas: P^ s oc /i 4 . The Wigner 
pressure is depicted in Fig. 5. 

With Eqs. (10), (16), (17) and (26), one has the zero- 
temperature EoS for the rainbow-truncation of the model 
defined by Eq.(4). The thermodynamical energy density 
follows: £t(/-0 = P" n(fJ-) — Pt(m)- NB. In assuming that 7] 



5 



is /z-independent, we neglect quark feedback on the gluon 
vacuum polarisation. In-vacuum, such effects are modest 
[28] . There is currently no reason to expect otherwise in- 
medium. We therefore anticipate that this feedback will 
have no qualitative impact on our results but may induce 
minor quantitative changes; e.g., a small reduction in fj, cr . 

The total pressure in the confining Nambu phase re- 
ceives a contribution from hadrons; viz . , P^ — ■ 
However, P^f is omitted by the rainbow truncation, 
which thus precludes an internally consistent descrip- 
tion of the transition between hadron and quark matter. 
Physically on the domain of Nambu-phase stability, one 
has dP^/dn > —dP^/d^; i.e., the vacuum rearrange- 
ment energy cost, evident in Fig. 2, is balanced (at least) 
by the gain from the response of hadron properties. This 
can happen whilst maintaining P^f /P^ <C 1. 

Nevertheless, for immediate application to compact 
astrophysical objects, our EoS must be augmented by 
a model for the nuclear matter EoS, which is the ac- 
tive branch for small chemical potential. Hitherto, a 
transition to quark matter was typically effected by 
a Maxwell construction and occurred when the pres- 
sures of the hadron and quark phases were equal; e.g., 



Refs. [29, 30, 31, 32]. However, our understanding of the 
nature of the vacuum pressure obviates the need for this 
prescription. Following a definition of the mapping be- 
tween quark and nucleon chemical potentials, we would 
execute a transition to deconfined and chirally symmet- 
ric quark matter at /z cr , with continuous pressure but 
discontinuous baryon number. 

One might ask whether our quark matter EoS should 
be preferred over those derived from bag- or NJL-likc 
models? The answer is contained in the following obser- 
vations. The model we have explicated is: (1) Poincare 
covariant, symmetry preserving, exhibits both confine- 
ment and DCSB, and provides a good description of in- 
vacuum hadron properties; and (2) exhibits coincident 
deconfinement and chiral symmetry restoring transitions 
at nonzero temperature and chemical potential. Neither 
(1) nor (2) can be said of bag- or NJL-like models. 
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